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ABSTRACT 

We describe a new approach to constrain the amplitude of the power spectrum of matter perturbations in the 
Universe, parametrized by cr% as a function of the matter density fio- We compare the galaxy cluster X-ray 
luminosity function of the REFLEX survey with the theoretical mass function of Jenkins et al. (2001), using the 
mass-luminosity relationship obtained from weak lensing data for a sample of galaxy clusters identified in Sloan 
Digital Sky Survey commissioning data and confirmed through cross-correlation with the ROSAT all-sky survey. 
We find cr 8 = 0.38 n>-°- 48+a27n °, which is significantly different from most previous results derived from comparable 
calculations that used the X-ray temperature function. We discuss possible sources of systematic error that may 
cause such a discrepancy, and in the process uncover a possible inconsistency between the REFLEX luminosity 
function and the relation between cluster X-ray luminosity and mass obtained by Reiprich & Btihringer (2001). 

Subject headings: cosmological parameters — galaxies: clusters — methods: analytical — surveys 



1. INTRODUCTION 

The present-day number density of galaxy clusters remains 
one of the most powerful constraints on the amplitude of matter 
perturbations in the Universe. This is usually defined in terms 
of as, which is the dispersion of the mass field smoothed on a 
scale of 8/T 1 Mpc -1 , where h is the present value of the Hubble 
parameter in units of 100 km s -1 Mpc -1 . Most often the local X- 
ray cluster temperature function is used for this purpose, given 
that the X-ray temperature has been the best observable from 
which to estimate cluster mass. While observation can readily 
give the X-ray temperature of clusters, theory can only easily 
predict the cluster mass function. To bridge the gap, theoretical 
modelling of clusters is used to provide a relation between mass 
and temperature, which in the most sophisticated treatments is 
taken to depend on both redshift and the underlying cosmology. 

However, a drawback in the use of the local X-ray cluster 
temperature function is that only a few tens of clusters have had 
their X-ray temperature estimated. Consequently, authors using 
samples that only partially overlap have obtained significantly 
different cluster temperature functions, and thus estimates for 
o-g (e.g. compare Eke, Cole & Frenk 1996; Henry 1997, 2000; 
Markevitch 1998; Blanchard et al. 2000; Pierpaoli et al. 2001). 
Statistically it would be more robust if the number density of 
local clusters could be estimated from much larger samples. 

We describe here a new approach which avoids working with 
the X-ray temperature function. We will instead use the lu- 
minosity function from the ROSAT-ESO Flux Limited X-Ray 
(REFLEX) galaxy cluster survey (Bohringer et al. 2001a,b), 
which contains 452 clusters, to estimate the local cluster num- 
ber density. We relate X-ray luminosity to cluster mass by 
taking advantage of a weak shear lensing analysis (Sheldon et 
al. 2001) of a sample of 42 galaxy clusters identified in data 
from the commissioning phase of the Sloan Digital Sky Sur- 
vey (SDSS), and cross-checked via correlation with the ROSAT 



All-Sky Survey (RASS) (Nichol et al. 2001). This relation is 
then used to compare the REFLEX luminosity function with 
the cluster mass function, leading to the estimation of a%, with- 
out any need to model the cluster mass-temperature relation (as 
in e.g. Viana & Liddle 1996, 1999). 

2. METHODOLOGY 

We consider spatially-flat cosmological models, where part 
of the energy density may be due to a cosmological constant, 
containing a spectrum of primordial adiabatic density pertur- 
bations. This family includes the current standard cosmologi- 
cal model, with a present-day matter density Q,q ~ 0.3, which 
provides the best fit to the full compilation of structure forma- 
tion data (e.g. Durrer & Novosyadlyj 2001; Wang, Tegmark & 
Zaldarriaga 2001). We take the present-day shape of the mat- 
ter power spectrum to be well approximated by that of a cold 
dark matter model with scale-invariant primordial density per- 
turbations and shape parameter V in the range [0.08,0.28]; this 
interval is the average of the best-fit values coming from the 
preliminary analysis of the 2dF (Percival et al. 2001) and SDSS 
(Dodelson et al. 2001; Szalay et al. 2001) data, taking into ac- 
count both statistical and systematic uncertainties. We will as- 
sume that T has an equal probability of taking any value within 
the interval given. Note, however, that ultimately the local clus- 
ter number density depends significantly only on cr% and 

The expected halo mass function for each set of cosmologi- 
cal parameters is estimated via the fitting function of Jenkins et 
al. (2001), obtained by analysing data assembled from various 
large A/-body simulations. Following White (2001) we con- 
sider the halo mass to be that given by the virial relation. The 
high-mass end of the halo mass function is better estimated us- 
ing the Jenkins et al. (2001) result rather than the usual Press- 
Schechter ansatz (Press & Schechter 1974). The latter leads to 
a smaller number of high-mass halos at fixed crs, thus requiring 
systematically higher values of a%, by about eight percent, in 
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order for the local cluster abundance to be reproduced. We cau- 
tiously model the uncertainty in the normalization of the mass 
function by means of a Gaussian distribution with a 10 percent 
dispersion around the mean (see Fig. 8 of Jenkins et al. 2001). 

The local cluster number density was obtained by integrat- 
ing the REFLEX X-ray luminosity function in the [0.1,2.4] keV 
band upwards from the mean luminosity of the 42 clusters in the 
SDSS/RASS sample, which is (0.17 ±0.03) x 10 44 /T 2 erg/s. 
Their average redshift is 0.1, which roughly coincides with the 
mean redshift of the REFLEX clusters with similar or higher lu- 
minosity. For the 42 SDSS/RASS clusters, Sheldon et al. (2001) 
obtained, via statistical weak lensing analysis, a mean projected 
mass within rgoo of (0.9 ± 0.2) x 10 14 h~ x M , where r 5 oo is the 
radius within which the cluster mean density falls to 500 times 
the critical density at the redshift of observation. Following 
Sheldon et al. (2001) we assume the cluster density profile to 
behave like that of a singular isothermal sphere, i.e. p(r) oc r~ 2 . 
We can now calculate the mean cluster mass within the three- 
dimensional radius rsoo, and then convert this to a mean virial 
mass. Taking this mass as the lower limit in the integration of 
the Jenkins et al. (2001) mass function will then yield the ex- 
pected local cluster number density for clusters more luminous 
than the mean of the SDSS/RASS sample, as a function of erg 
and Clo- By comparison with the REFLEX estimate, best-fit 
values for erg as a function of flo can then be obtained. 

3. RESULTS 

Table 1 shows the 95 per cent confidence interval on erg ob- 
tained using the REFLEX X-ray luminosity function and the 
relation between cluster mass and X-ray luminosity for the full 
SDSS/RASS sample. This interval was determined via Monte 
Carlo simulations, which incorporated all the uncertainties pre- 
viously mentioned that affect the present estimation of erg. We 
find that the most probable value for erg can be accurately rep- 
resented by the fitting function 

<jz = 038n- QA ^ y2m \ (1) 
with a 95 per cent uncertainty around 15%. 

This result is significantly lower than, and barely compati- 
ble with, most determinations of erg based on the local clus- 
ter X-ray temperature function (e.g. compare with Eke, Cole & 
Frenk 1996; Henry 1997, 2000; Viana & Liddle 1999; Blan- 
chard et al. 2000; Pierpaoli et al. 2001), with the exception of 
Seljak (2001). His analysis differs from the others in that he 
used the relation between cluster temperature and mass derived 
by Finoguenov, Reiprich & Btihringer (2001) from X-ray data, 
rather than one obtained from hydrodynamical A^-body simula- 
tions. The earlier work of Markevitch (1998) had already hinted 
at lower values for erg, at least in the case of Slo — 0.3, if actual 
X-ray data was used to relate cluster temperature to mass. Sim- 
ilar results have been reached by Borgani et al. (2001), based 
on ROSAT Deep Cluster Survey (RDCS) data and the observed 
cluster X-ray temperature to luminosity relation, and Reiprich 
& Btihringer (2001), by means of an empirical cluster mass 
function derived using X-ray data from a large cluster sample. 

Other methods of measuring erg lead to conflicting results. 
While high-redshift Lyman-a forest analyses (Croft et al. 1999, 
2000; McDonald et al. 2000), support our findings, estimates 
based on cosmic shear data (Htikstra, Yee & Gladders 2001; 
Maoli et al. 2001; Van Waerbeke et al. 2001) tend to favour 
higher values for erg than those obtained here. Also, the 2dF 
galaxy survey, when combined with measurements of the am- 
plitude of temperature anisotropics in the cosmic microwave 



Table 1 

Matter power spectrum normalization using the 
reflex luminosity function and the full sdss/rass 
sample, shown with 95% confidence error bars. 
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background radiation (Efstathiou et al. 2001 ; Lahav et al. 2001), 
seems to require a value for erg close to what we have found. 

In order to test whether the low values obtained for erg could 
be due to hidden systematic errors in the weak lensing method 
used for cluster mass estimation, we repeated the calculation 
of erg for the two SDSS/RASS sub-samples discussed within 
Sheldon et al. (2001), namely 

1 . The 27 clusters with the lowest X-ray luminosities, on 
average (0.09 ±0.02) x 10 44 /T 2 erg/s. They have a mean 
redshift of 0.09 and a mean projected mass within rgoo 
ofM 50 o(l) = (0.7±0.2)xl0 14 r 1 M Q . 

2. The 15 clusters with the highest X-ray luminosities, on 
average (0.51 ±0.04) x 10 44 /T 2 erg/s. They have a mean 
redshift of 0.17 and a mean projected mass within r$oo 
of M 50 o(2) = 2.7+^ x 10 14 h~ l M . 

Given that most of the REFLEX clusters with similar luminosi- 
ties to those in the SDSS/RASS sample have a redshift between 
0.05 and 0.2 (Btihringer et al. 2001a), with the higher luminos- 
ity clusters typically being at higher redshifts, we will assume 
that the REFLEX luminosity function provides a good repre- 
sentation of the underlying cluster luminosity function over this 
redshift interval. This is supported by an analysis of the Bright- 
est Cluster Survey (BCS), which showed that there is no strong 
evidence for evolution in the cluster luminosity function out to 
at least z = 0.2 (Ebeling et al. 1997). 

Surprisingly, the most probable values for erg according to 
each of these two sub-samples are rather different. Sub-sample 1 
yields erg = 0.37 C2- a52+0 - 32O » ? almost indistinguishable from the 
result obtained from the full SDSS/RASS sample, while sub- 
sample 2 gives erg = 0.50 ft" ' 47 * ' 24 ^ which is substantially 
higher. This latter result is much more in line with standard 
erg estimates from the local cluster X-ray temperature function. 

A. DISCUSSION 

In this section, we discuss our results and present several tests 
to assess their robustness. First, to determine whether our re- 
sults indicate an internal inconsistency in the weak lensing anal- 
ysis of Sheldon et al. (2001), we performed Monte Carlo simu- 
lations where we used the data from sub-sample 1 to compute 
the inferred mean projected masses within f-500 for the clusters 
in sub-sample 2 and vice-versa, for various values of We in- 
cluded all the uncertainties previously mentioned in the simula- 
tions. We found that, provided f2o > 0.3, sub-sample 1 implies 
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a value for M5oo(2) at least as large as that deduced by Sheldon 
et al. (2001) less than 5 per cent of the time, while sub-sample 
2 implies a value forM5oo(l) as small as that deduced by Shel- 
don et al. (2001) less than 10 per cent of the time. For values of 
Slo between 0.1 and 0.3, the discrepancy decreases with Slo, but 
not significantly. We therefore conclude that, within the con- 
text of the cosmological models we discuss, the mean projected 
masses for the two SDSS/RASS cluster sub-samples presented 
in Sheldon et al. (2001) are only barely compatible within the 
uncertainties associated with the estimation of ag. 

This inconsistency could result from the assumption that the 
cluster mass profile is that of a singular isothermal sphere (SIS), 
which was motivated by the fact that within the radius for which 
there was weak lensing data, the shear profile was close to that 
of a projected SIS. Although observationally the issue is still 
unresolved (e.g. Irwin, Bregman & Evrard 1999; Irwin & Breg- 
man 2000; White 2000; De Grandi & Molendi 2001; Komatsu 
& Seljak 2001), it has become clear that A^-body simulations 
produce clusters which on average have outer mass profiles, 
up to the cluster virial radius, that behave more like p{r) oc 
r" 2 5 (Navarro, Frenk & White 1995, 1996, 1997; Thomas et 
al. 1998, 2001; Tittley & Couchman 2001). Consequently, we 
investigated the effect on the virial mass to X-ray luminosity 
relation of assuming an outer cluster mass profile parametrized 
by p(r) oc r _/3 . We varied (5 between 2 (SIS case) and 3, and, 
as before, we took the mean projected masses within r^Q to 
be those given by Sheldon et al. (2001). This translates to 
fixing the inner projected cluster masses to be M5oo(l) within 
0.39 /T 1 Mpc for sub-sample 1, M 50O (2) within 0.57 IT 1 Mpc 
for sub-sample 2, and 0.9 ± 0.2 x 10 14 M Q within 0.41 /T 1 Mpc 
for the full SDSS/RASS sample (Sheldon et al. 2001). As ex- 
pected, we find that as (3 increases, the cluster virial mass de- 
creases, hence lowering a%. This change is greater for smaller 
Slo- Consequently, we find that the value for erg derived using 
the Sheldon et al. (2001) and REFLEX data is barely affected 
by assuming different outer mass profiles if Slo is close to one, 
whereas for low values of Slo assuming a outer mass profile dif- 
ferent from the true one may introduce a significant systematic 
error in the calculation of a%. For example, for f3 = 3 we obtain 
erg = 0.37Slo°- 41+0 ' 240 ° as the most probable value that results 
from the weak lensing analysis for the full SDSS/RASS sam- 
ple. However, it turns out that changing the outer cluster mass 
density profile does not significantly mitigate the discrepancy 
between the values obtained for erg using the two SDSS/RASS 
sub-samples, even for low values of Slo- 

The discrepancy between the erg results obtained for the two 
Sheldon et al. (2001) sub-samples could be alleviated if the 
mean projected mass of the high-luminosity sample is overes- 
timated, or that of the low-luminosity sample underestimated. 
The first hypothesis is much more probable, and could be due to 
a contribution to the mean projected mass by filamentary mate- 
rial infalling into the clusters along the line of sight (Cen 1997; 
Metzler et al. 1999; Reblinsky & Bartelmann 1999; Metzler, 
White & Loken 2001). Alternatively, the mean X-ray lumi- 
nosity of the high-luminosity sample may have been underesti- 
mated, or that of the low-luminosity sample overestimated. The 
latter could be caused by AGN contamination, but the initial 
study of Miller et al. (2002, in preparation) indicates that this 
is not a serious problem in the case of the SDSS/RASS cluster 
survey. Furthermore, the observed relationship between X-ray 
cluster luminosity and measured velocity dispersion (from the 
SDSS spectroscopic sample) for both the Sheldon et al. (2001) 



and the Miller et al. (2002, in prep.) samples are in excellent 
agreement with the Lx-cr v relation of Mahdavi & Geller (2001). 

Reiprich & Bohringer (2001) used ROSAT and ASCA X-ray 
data on 106 clusters to obtain the relation between X-ray lumi- 
nosity in the [0.1,2.4] keV band and cluster mass in the form of 
M 5 oo, with the cluster masses estimated assuming hydrostatic 
equilibrium (Finoguenov et al. 2001; Reiprich & Bohringer 
2001). They found L x = 10" 18 ' 59±1 - 23 x M^ 75±0 084 where L x 
is in units of 10 40 erg/s and M500 in units of M Q (assuming h = 
0.5). Substituting the mean luminosities for the SDSS/RASS 
full sample and two sub-samples into this relation, one finds 
that the Sheldon et al. (2001) estimates forMsoo are well within 
the (extremely wide) range of possible values. Conversely, the 
relation between Lx and M500 can be estimated by combining 
the Sheldon et al. (2001) results with the shape of the cluster 
mass function, assumed well described by that of Jenkins et 
al. (2001), and the data on the cluster luminosity function from 
the REFLEX survey (Bohringer et al. 2001a,b). Assuming the 
Lx -M500 relation to be a power-law, we found that the normal- 
ization is essentially defined by the Sheldon et al. (2001) data, 
as expected, while the exponent is mainly governed by the rel- 
ative shape of the mass and luminosity (cumulative) functions. 
Taking into account all the uncertainties involved in the normal- 
ization of the Jenkins et al. (2001) mass function by means of 
the Sheldon et al. (2001) data for the full SDSS/RASS sample, 
and those associated with the REFLEX luminosity function, we 
obtained through Monte Carlo simulations a 95 per cent confi- 
dence interval of [1.6,3.7] for the exponent of the Lx-Msm re- 
lation. The allowed interval does not change significantly if 
the data for either of the two SDSS/RASS sub-samples is used 
instead to normalize the mass function. The low-luminosity 
sample yields [1.8,3.2], while from the high-luminosity sam- 
ple we get [1.5,3.1]. In any case, surprisingly, the preferred 
value is close to 2.2, substantially higher than the 1 .58 ± 0.08 at 
lcr found by Reiprich & Bohringer (2001). This analysis was 
performed for a flat universe with Slo = 0.3, though very similar 
results were found for Slo = 1 . In order to confirm that the expo- 
nent of the L x -M500 relation is only weakly determined by the 
normalization of the mass function, and thus by the Sheldon et 
al. (2001) data, we varied erg between 0.5 and 1.2 (for Slo = 0.3) 
and found that the preferred value for the exponent changes 
from 2.4 to 1.7. Therefore, assuming the Jenkins et al. (2001) 
mass function provides an accurate description of the cluster 
mass function, the discrepancy just found on the exponent of 
the Lx -M500 relation means that such a relation as obtained in 
Reiprich & Bohringer (2001) is at best only marginally consis- 
tent with the Sheldon et al. (2001) data taken together with the 
REFLEX luminosity function. And only if erg is at the higher 
end of recent estimates (see e.g. Viana & Liddle 1999) can the 
Reiprich & Bohringer (2001) L x -Msqo relation be made con- 
sistent with the REFLEX luminosity function, within the con- 
text of the cosmological models discussed in this paper. 

In calculating erg we assumed that, on average, a galaxy clus- 
ter with the mean X-ray luminosity of the SDSS/RASS full 
sample has M500 equal to the value estimated in Sheldon et 
al. (2001) for such a sample. This assumption is prone to sev- 
eral biases, one being that, assuming that there is some disper- 
sion in the associated luminosity, the clusters of a given mass 
that preferentially end up in a sample selected in the manner of 
the SDSS/RASS are the most luminous ones. This bias leads 
to an underestimation of the correct mass corresponding to a 
given luminosity, and consequently to an underestimation of erg, 



4 



P. T. P. Viana, R. C. Nichol and A. R. Liddle 



though it is difficult to say by how much. The same effect takes 
place due to another, more subtle, type of bias, arising from 
the fact that the L x -Msqo relation is not linear. As a result, for 
a given sample the mean M500 is not proportional to the mean 
Lx- Given that M500 oc L\ with n < 1, assuming that a cluster 
with the mean luminosity has a mass M500 equal to the sample 
mean underestimates M500 for such a cluster (because actually 
the mean M500 is proportional to the mean of L x ). We have 
found that, for the luminosity dispersion of the clusters in the 
SDSS/RASS full sample, the most probable value of M500 for a 
cluster with the mean L x may be underestimated by around 20 
per cent. This percentage is robust to changes in the assumed 
value for the exponent n (between 0.3 and 0.7) and to the pos- 
sibility of dispersion in the Lx-Msm relation. This bias leads 
to a possible underestimation of erg close to 10 per cent, almost 
independent of However, given that the mean M500 for the 
SDSS/RASS sample is in fact not a mean of several indepen- 
dently calculated M500, but the result of a mean shear profile, 
and that the mean Lx is weighted by each cluster contribution 
to that profile, the above considerations may not be directly ap- 
plicable to the case in hand. Note that because the dispersion in 
luminosities is smaller for either of the two SDSS/RASS sub- 
samples as compared to the full sample, the underestimation in 
erg is smaller when it is estimated from the sub-samples, being 
closer to 5 per cent. Unfortunately, none of these biases seem 
to be able to significantly narrow the discrepancy between the 
erg values derived using each SDSS/RASS sub-sample. 

5. CONCLUSIONS 

We have applied a new approach for constraining the normal- 
ization of the matter power spectrum, using the REFLEX X-ray 
luminosity function and the relation between cluster X-ray lu- 
minosity and mass obtained through weak lensing data from 
a preliminary small sample of SDSS/RASS clusters. We ob- 
tained erg values significantly lower than other estimates based 
on cluster abundance data, with the exception of the recent re- 
sults by Borgani et al. (2001), Reiprich & Bohringer (2001) 
and Seljak (2001). However, systematic biases may affect our 
analysis, given that barely consistent results are obtained when 
using subsets of the weak lensing data. This may be due to the 
small sample of clusters used, or an artifact of the techniques 
used in Sheldon et al. (2001) for co-adding clusters to produce 
an ensemble averaged weak lensing signal. In the process, we 
found that comparing the REFLEX luminosity function and the 
Jenkins et al. (2001) mass function implies that the relation be- 
tween X-ray luminosity and cluster mass may be significantly 
steeper than previously thought. 

The SDSS/RASS data set we have used will be dwarfed by 
the final SDSS/RASS catalogue (see Nichol et al. 2001), and 
surveys using the XMM-Newton satellite should supply much 
greater information on cluster luminosities (e.g. Romer et al. 
2001). The prospect of considerably improving the constraint 
on erg using this approach in the future is therefore great. 
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